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Introduction 

o 



The definition of a partially invariant solution (PIS) as a natural generalization of an 
invariant solution (IS) of differential equation was first suggested by L.V. Ovsiannikov 
[21 II]- At present there are known many examples of partially invariant solutions, 
mostly for models of fluid mechanics [3]-[T7]. In contrast to the invariant solution, the 
partially invariant one is described by an overdetermined system of equations. The 
main difficulty of PIS construction is the investigation of compatibility and involutivity 
[18] of this system. This complicated analysis requires individual approach for every 
particular solution. On the other hand, partially invariant solutions usually possess 
higher arbitrariness then the invariant ones, and, therefore, are more interesting for 
physical interpretation. 

In the present work it is noticed that the majority of PISs of a given system of 
equations can be obtained as a result of the two-steps construction: As the invariant 
reduction of some partially invariant solution of the higher rank. This introduce a 
hierarchic structure in the set of all partially invariant solutions of a given system 
of equations. Those PISs which do not have a nontrivial two-steps representation 
are referred to as indecomposable ones. Analysis of all indecomposable PISs allows 
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consequent construction of the whole set of partially invariant solutions by means of 
invariant reductions only The latter do not require investigation of overdetermined 
systems of equation and, thus, are simpler. This investigation supplemented with LOT 
lemma FTJ] introduces the hierarchic structure on the whole set of the group-invariant 
submodels of a given system of differential equations. 

The theory of partially invariant solutions taken from [1] is briefly recounted in 
sections 1 and 2. Hierarchy of partially invariant submodels is discussed in sections 3 
and 4. Example of PISs hierarchy for shallow water equations is given in section 5. 
Section 6 is devoted to the complete classification of indecomposable regular PISs for 
ideal MHD equations. 

1. Partially invariant manifolds 

Let G r = {T a : x = f(x,a)}, a G A C W be a local Lie group of transformations 
acting in space x G MJ 1 . Basis of the corresponding Lie algebra L r of infinitesimal 
generators is given by X a = ^(x)d x i, (i — 1 ; n; a — 1, r). Hereafter the Einstein 
summation convention on the repeating upper and lower indices is adopted. Let us 
observe a manifold M. regularly defined by equations 

M : 4> a (x) = Q, a = l,...,s; rank || dip a /dx k \ | = s. (1.1) 

Hereafter rankM(x) denotes the maximal rank of matrix M(x) for various values of x. 

Definition 1 Orbit of point x under G r group action is a set of points 0(x) = 
{f(x,a) | a G A C W}. Orbit 0(M) of the manifold M. is the locus of orbits of 
all points x G A4: 

0(M) = {f(x, a) | x G M, a G A C R r } 

Let us introduce the following integer characteristic. 

Definition 2 Defect 5 (M,G r ) of the manifold A4 under G r group action is a difference 
between the dimensions of the orbit 0(M) and of the manifold M. itself: 

5{M,G r ) = d\mO(M)-dimM. (1.2) 

Defect of the manifold is an important characteristics showing the degree of non- 
invariancy of the manifold M. under the action of G r . Defect of an invariant manifold 
is zero because orbit of each point of the manifold lies entirely on the manifold. 

Orbit of an arbitrary manifold under the group action is itself an invariant manifold 
of the group because, by definition, it contains orbits of all its points. Moreover, orbit 
0(M) is the minimal invariant manifold of the group G r , containing Ai. Thus, it can 
be described in terms of the functional relations between the invariants of the group. 
Let the complete set of functionally independent invariants of G r be chosen in the form 
/ = (^(x), . . . , I t (x)), where t = n — r*, and r* = rank ||£„(x)||. By virtue of the 
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theorem of representation of a non-singular invariant manifold pQ the orbit of Ai can 
be written in the form 

^(l\x),...,I t (x)) = 0, t = 1,...,L (1.3) 

Definition 3 Under the condition of regularity of specification / TOj) . i.e. rank ||<9$ T j dl \ 
I the number 

p(M,G r )=t-l (1.4) 

is referred to as the rank of the partially invariant manifold M. with respect to the group 
G r . Pair of integers (p, 5) define the type of the partially invariant manifold M.. 

Rank of the manifold is equal to the dimension of the orbit O(Ai) in space of invariants 
of the group G r . In practical calculations formula (11.41) is inconvenient, because it rely 
on the invariant representation ( 11. 31) of the orbit of the manifold At, which may be not 
known explicitly. However, by using ( II. 2p . the rank can be found in terms of codimension 
s of the initial manifold Ai and of its defect 5(A4, G r ): 

p(M,G r )=5(M,G r )+t-s. (1.5) 



2. Partially invariant solutions 

Let us observe a system of differential equations 

E : F a {x,u,u, ...,u) = 0, a=l,...,s. (2.1) 

1 k 

The main space is Z = M. n (x) x M m (w). By u we denote the set of all p-ih order 

v 

derivatives: | dx^dx^ } ' ^y s ^ em P-ip admit Lie group G r = {T a : Z x W — > Z}. 
Action of G r can be prolonged on the derivatives in usual manner [H [20]. Let Lie 
algebra of infinitesimal generators of G r be 

L r = {X a = £ l a (x,u)d x i +r]*(x,u)d u k, a = l, 

For a fc-dimensional subalgebra H C L r matrix 

/ £(x,u) ... g(x,u) \ 



,r}. 



H(0 



and extended matrix 



( £{(x,u) ... ii{x, u) 7]\(x,u) ... rf^{x, u) \ 



\£l(x,u) ... £,%(x,u) r)l(x,u) ... rf£{x,u) ) 



are introduced. Let 

U : u l = ip l (x), i = l,...,m. 
is a solution of equations (12. ip . 



(2.2) 
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Definition 4 Solution U is referred to as H -invariant solution (H-IS) of the system of 
equations E if the manifold U C Z is an invariant manifold under the subgroup H C G r 
action. 

The necessary condition of H-IS existence is convenient to formulate in terms of the 
corresponding Lie subalgebra of infinetesimal generators. 

Lemma 1 Lie subalgebra H C L r generate H -invariant solution of the system E if the 
following equality holds 

rank H{£) = rank H(£, rj). (2.3) 

Generalization of the notion of the invariant solution leads to the following. 

Definition 5 Solution U is called H -partially invariant solution (H-PIS) of the system 
E if the manifold U C Z is a partially invariant manifold under H C G r action. 

PISs are usually constructed on subalgebras of the admissible algebra which do not 
satisfy the necessary condition (12.31) . Definitions of rank and defect of a partially 
invariant manifold are transferred naturally on PISs. However, there is some specifics 
owing to the distinction of roles of x and u variables. 

Let us introduce the following integer characteristics of H: 

t = m + n — rankiJ(£, rj) — the total number of invariants of H; 
a = n — rank H (£) number of invariants of group H, which depend 

only on x; 

ji — t — a — number of invariants essentially depending on u. 

Partially invariant solution specified in form of a manifold $ by formulae (12.21) 
have codimension s = m. Rank of manifold $ can be calculated by formula (jl.5p as 
p = 5 + t — m. Orbit of manifold $ is an invariant manifold of group H, therefore 
it may be specified by some functional relations of the form (I1.3P for the invariants of 
the group. For the sake of simplicity we assume that the invariants are chosen in the 
separated form 

f I 1 {x,u),...,I' 1 {x,u), 
I : <^ 2.4 
[ A 1 = I t -^ 1 {x),...,X a = I t {x). 

At that 

a(jH s «),...,i»(,,«)) (2 . 5) 

d^ 1 , . . . ,u m ) 

Let us form equations of the orbit of a partially invariant manifold $ under the 
group H action. At that it is required to specify the rank of a partially invariant 
manifold. This can be taken as any integer p, satisfying the inequality 

a < p < min{n, t}. (2-6) 
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Definition 6 Partially invariant solution is called regular if p — a. 

Equations of the orbit 0(<&, H) are constructed as the following set of functional relations 
between the invariants (j2.4p of the group H 

$ r (A 1 ,...,A ff , = 0, t = l,...,t- p; rank || d^/dl k \\ = t-p (2.7) 

with unknown functions $ T . In practical calculations this relations can be taken in the 
resolved form, e.g. 

F = ^(A\ ...,\°, ...,/'*), r=l,...,t-p. (2.8) 

Although, this form is more convenient for PIS computation, it is not unique in case of 
irregular solutions (i.e., p > a). Therefore, in the theoretical analysis it is preferable to 
refer to the general form ( 12. 71) . 

The dimension of the manifold ( 12.71) in the space of invariants is p. Equations 
(12. 7p can be solved with respect to t — p functions u by virtue of conditions (12. 5p . 
(12.61) . Without loss of generality one can assume these functions to be w 1 , ... , u t ~ p . The 
remaining 5 = m — t + p functions m'~ p+1 , . . . , u m do not have representation in terms of 
invariants of H and initially are not restricted by any extra assumptions. Thus, there 
appear 5 non-invariant functions, which are assumed to depend arbitrarily on x. 

u l - p+1 =w 1 {x),...,u m = w 5 {x). (2.9) 

Formulas ( 12. 7\\ . ( 12.91) define the representation of the partially invariant solution of the 

type (p,S). 

Remark 1 In practical calculations it is necessary to observe all n on- equivalent 
possibilities for solution of the orbit equations {2.1 ) with respect to functions u and 



obtaining of the representation Ii2.9\) for non-invariant functions. In what follows the 



representation of solution will refer to the combination of the orbit equations \2. 1\ ) with 
all possible representation of non-invariant functions of the form Ii2.9\) . 

Substitution of the representation ( 12.71) . (12.91) of solution into the system of 
equations (12. ip leads to a factor-system of differential equations for the invariant 
functions <3> fc , k = 1, . . . , p, and non-invariant functions wi , j = 1, . . . , 5. The factor- 
system of a partially invariant solution contains a subsystem E/H for invariant functions 
and invariant variables, and equations II for the non-invariant functions. System II 
of equations should be observed as an overdetermined system for the non-invariant 
functions wK At that, all the invariant functions $ fc are assumed to be known 
from a solution of the invariant subsystem E/H. The compatibility conditions of II 
usually extend both the invariant part E/H, and the system II itself. The purpose 
of investigation at this stage is to bring system IT to involution, i.e. to obtain all its 
compatibility conditions to prove its self-consistency. Unfortunately, it is impossible 
to trace this process in general form. If this step is performed, the factor-system 
finally takes the form of a union of a subsystem E/H for the invariant functions and 
of compatible on the solutions of E/H system II for determination of non- invariant 
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functions. This reduced system is simpler then the original system E because E/H 
involves less independent variables, and II contains less unknown functions. 

Definition 7 The union of the factor-system E/H and of the system U will be referred 
to as H -partially invariant submodel of the system of differential equations E. 

3. Partially invariant solutions hierarchy 

Let H,N C G r be subgroups, such that H is a normal divisor in N: H < N. Suppose 
that H does not satisfy conditions (12.31) of invariant solution existence. Let if -partially 
invariant submodel of E be known. In the sequel the following question is investigated: 
under which conditions on group N there exists N/ H-invariant solution for H-PIS and 
how does it relate to iV-PIS of equations El 

Owing to the one-to-one correspondence between local Lie groups of transforma- 
tions and their Lie algebras of infinitesimal generators [H ED], later on we do not dis- 
tinguish between these two objects, denoting them by the same letter. All facts below 
proven in the Lie group language can be translated into the Lie algebras language and 
vice versa. 

Lemma 2 The factor group NJ H have induced action in the space of invariants of the 
group H . 

Proof. Factor group N/H is a set of all left equivalence classes gi=goH = {goh\hE 
H}. Let J be an invariant of the group H action. Action of gi on J is defined 
as gi(J) —go h(J) = g(J), and, obviously, does not depend on the choice of the 
representative h of the equivalence class. According to the condition H < N one has 
g~ 1 Hg C H for each g G N. Let us show that if J is invariant with respect to H 
then with every g G N function g(J) is also invariant under group H action. Indeed, 
according to g^ 1 oho g( J) = hi(J) = J we have h(g(J)) = g(J) for each h G H. Thus, 
the factor group N/H action is defined and closed in the space of invariants of group 
H. □ 

Lemma 3 Group N/H is admitted by the system of differential equations E/H; group 
N is admitted by the system U on solutions of E/H. 

Proof. Let us show that normal extension of group H acts on the set of ii-invariant 
manifolds. Indeed, let H < N and Ai is some if-invariant manifold. For any h G H and 
g G N there exists hi G H such that h o g = g o h\. Hence, 

h(g(M)) = g(h 1 (M))=g(M). 

Thus, the transformed manifold g{M) is also if- invariant. 

Orbit of H-PIS, given by equations (12. 7ft . is an invariant manifold of the group H. 
Then, any transformation of group iV transforms an orbit of H-PIS into an orbit of H- 
PIS. The orbit of an arbitrary H-PIS is determined by the system E/H. Action of the 
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factor-group N/ H is closed in the space of invariants of group H . Hence, transformations 
of factor-group N/H acts on the set of solutions of E/H system, i.e. are admitted by 
this system of differential equations. Next, group iV acts on the set of solutions of the 
original system N . Besides, it conserves an orbit of a solution in class of H-PIS. Thus, 
N acts in the set of if-PISs, i.e. is admitted by system II on solutions of system E/H. 
□ 

The necessary condition of iV/f/- invariant solution existence is convenient to 
formulate in the Lie algebraic language using infinitesimal generators of the observed 
Lie group. 

Lemma 4 For N / H -invariant solution of the factor-system E/H to exist the following 
condition should be satisfied 

rankiV(£,77) - rankiV(f) = rank H (£,77) - rank#(£). (3.1) 

Proof. Let us transform infinitesimal generators of H into the coordinate system, which 
flattens the group action: 

y 1 = A 1 (a;), . . . , y° = A CT (x); y° +l = x° +1 , ...,y v " 

v 1 = I x (x, u), . . . , f M = /^(x, u); v^ +1 = it M+1 , . . . , v r> 

Without loss of generality, this transformation is non-degenerate. Group H acts 
transitively in the space of variables (y a+1 , ■ ■ ■ , y n ] v . . . , v m ). Representatives of 
the factor-algebra N/H have the form X = X + Y, where Y G H. Operator Y does 
not contain any differentiations with respect to invariant variables y 1 , . . . , y a , v 1 , . . . , v^. 
Operator X by virtue of lemma [2] projects into the space of invariants of group H, i.e. 
its coefficients at the differentiations with respect to invariant variables do not contain 
non-invariant variables. Thus, operator Y does not participate in projection of operator 
X into the space of invariants of group H at all. The Lie algebra of projections {X} 
corresponds to the induced action of the factor-group N/ H in space of invariants of H. 
This Lie algebra should satisfy necessary conditions of the invariant solution existence 
( 12. 3p . In what follows we check this conditions explicitly in terms of coordinates of 
infinitesimal operators of Lie group N. 

Matrix of coordinates of infinitesimal generators of Lie algebra N in (y,v) 
coordinates have the block structure: 



x 



u 



(3.2) 





( 


A 





B 


\ 













C 






\ K 


L 


R 


s 


J 



First two columns of blocks (to the left of the vertical double line) correspond to 
differentiations with respect to independent variables y, the remaining two columns (to 
the right of the double line) are coordinates at differentiations with respect to dependent 
variables v. Number of columns in each block separated by vertical lines is equal to a, 
n — a, fi, and m — /j correspondingly. The matrix horizontal division is such that first 
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dim if rows (above the horizontal line) contain basic operators of Lie algebra H. The 
remaining operators (below the horizontal line) complete H to N. 

Blocks of matrix iV(£, 77) with coordinates (row, column) = (1, 1), (2, 1), (1,3), and 
(2, 3) are equal to zero because coordinates at differentiations with respect to invariant 
variables in infinitesimal generators of H vanish. Algebra H does not satisfy the 
necessary conditions of existence of an invariant solution (12.31) . Hence, by performing a 
suitable combination of rows above the horizontal line in matrix iV(£, 77), on can make the 
block with coordinates (2, 2) to be zero with rank A = n—a. As long as only the ranks of 
blocks are of interest, the linear combinations may be taken with coefficients depending 
on all variables. The resulting block C is also non-degenerate: rankC = m — fx. Thus, 
blocks A and C are square non-degenerate matrices of dimensions (n — a) x (n — a) and 
(m — //) x (m — //). 

Let us turn to the rows of N(£, rj) placed below the horizontal line and corresponding 
to basic elements of complement of H to N. By a non-degenerate combination of these 
rows with upper n — a rows of matrix iV(£, rj), one can zero block L. Next, by virtue 
of the non- degeneracy of block C , one can zero block S. At that, blocks K and R 
stay unchanged. Differential operators with coordinates from blocks K and R acting in 
the space of invariants of group H, are the required infinitesimal operators of induced 
action of the factor-group N/ H in the space of invariants. The necessary condition of 
existence of invariant N/H -solution (12.31) is that the rank of block K is equal to the 
rank of matrix composed from blocks K and R: 

rank isT = rank (K, R). (3.3) 

This condition can be reformulated in terms of ranks of blocks of the complete matrix 
N(£,rj). Let columns of blocks of matrix N(£,rj) be denoted by Roman numbers i - 
IV. Notice, that by virtue of non-degeneracy of blocks A and C columns II and TV 
are linearly independent with each other and with columns i and III. 

rank (i, II) = rank K + n — a, 

rank (i, II, III, IV) = rank (K, R) + n - a + m - fi. 
Subtraction of the first equality from the second one gives 

rank (I, II, III, IV) — rank (I, II) = rank (K, R) — rank K + m — fx. 
Relation (13.31) is satisfied if and only if 

rank (I, II, III, IV) - rank (i, II) = m- (i, 

which is equivalent to ( 13.11) condition. As long as change of variables ( 13. 21) does not 
affect ranks of matrices, this condition can be verified in initial coordinate system. □ 

Remark 2 The necessary conditions formulated in lemmas U\ and [7] guarantee the 
possibility of construction of representations of the corresponding invariant solutions 
and finding factor-systems of equations for invariant functions. These conditions are 
not sufficient because compatibility of the obtained factor-systems can not be proved a 
priory. 
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4. Two-step construction of the solution 

The information obtained about partially invariant solutions allows formulating the 
following statement. 

Theorem 1 Let system of differential equations E admit Lie group of continuous 
transformations N. Suppose that there is a normal divisor H in N , which does 
not satisfy the condition of existence of the invariant solution 112.5]) . but fulfils the 
requirement $3.1\) . Then, for the factor-system E/H, corresponding to H-PIS there 
exists an invariant solution with respect to the factor- group N/H. Moreover, the factor- 
system of the invariant solution (E/H) /(N/H) is equivalent to the factor-system of the 
partially invariant solution E/N. 

Proof. The possibility of construction of iV/if-invariant solution of if -PIS is already 
shown in lemmas EHU The only thing left to demonstrate is the equivalence of factor- 
systems E/N obtained directly and by using the two-steps method as (E/H) /(N/H). 

Suppose, that the dimension of the factor-group N/H is equal to k. Owing to the 
condition (I3.3P one can assume invariants A 1 , . . . , A CT , and I 1 , . . . , 7 M of the group H to 
be chosen in such a way that functions 

A 1 ,..., A'"", J 1 ,...,/" (4.1) 

form the basis of functionally independent invariants of the group N. This implies, that 
the action of the factor-group N/H is transitive in the subspace M K (A <T_K+1 , . . . , A CT ), 
whereas variables (14.11) are independent invariants of N/ H . The representation of N/ H- 
invariant solution of the factor-system E/H is obtained by the demand of independency 
of functions $ in ( 12.7ft on variables A°"~ K+1 , . . . , A CT . Exactly the same representation of 
the orbit of iV-PIS of equations E is obtained directly. Thus, representations of iV-PIS 
and of (E/H ) /( N/H)-1S coincide, hence the factor-systems followed are equivalent as 
well. □ 

Thus, in the set of partially invariant solutions of investigated system there is a 
hierarchic structure. 

Definition 8 Partially invariant submodel is called indecomposable if it can not be 
represented as the non-trivial combination of a partially invariant and invariant 
submodels. 

Investigation of only indecomposable submodels allows significant reduction of efforts in 
enumeration of all partially invariant solutions of a given system of differential equations. 
Indeed, the most labor-intensive step of involutivity analysis for overdetermined systems 
of differential equations should be done only for indecomposable submodels. The 
remaining submodels are obtained from the indecomposable ones by means of only 
invariant reductions, which is usually much simpler. The following two sections 
demonstrate this approach on examples of the shallow water and ideal MHD equations. 
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5. Shallow water equations 

The equations, describing motions of a thin water layer over a flat bottom are observed: 

u t + uu x + vuy + h x = 0, 

E : v t + uv x + vv y + h y = 0, (5.1) 

fk + {uh) x + (vh) y = 0. 

Here (u, v) is a particle's velocity vector, h is the depth of the water layer. The basic 
space here is ~R 3 (t,x,y) x M 3 (u, v , h), hence n = m = 3. The admissible algebra L 9 [2] 
is generated by operators (notations of paper [21] are adopted): 

X\ = d x , X 2 = d y , X 4 = td x + d u , X 5 = td y + d v , 
X 9 = xd y - yd x + ud v - vd u , X w = d u 
X n = xd x + yd y + ud u + vd v + 2hd h , 

X\2 = t 2 d t + txd x + tydy + (x — tu)d u + (y — tv)d v — 2thdh- 
X 13 = 2td t + xd x + ydy - ud u - vd v - 2hd h . 
Let us observe a partially invariant solution given by Lie subalgebra 
N = {X U X 4 , X 10 + X 12 }CL 9 . 

The subalgebra H = {X%, X 4 } in is selected. It is easy to check that H is ideal in 
N. Condition (13. ip is satisfied, therefore one can apply the two-steps algorithm. 
The complete set of functionally independent invariants of H is 

t, y, v, h. 

Here n = 3, m = 3, t = 4, a = 2, /i = 2. Let us construct an indecomposable if-PIS 
of rank 2. Equation of orbit of a partially invariant solution (12.71) can be written in an 
explicit form 

v = v(t,y), h = h(t,y). 

Defect of the solution is S = 1. There is only one non- invariant function u, which is 
supposed to depend on all independent variables: 

u = u(t,x,y). 

Substitution of the obtained representation of solution into the initial system (15.11) 
gives the submodel equations. The first and the third equations of system (15.11) form 
overdetermined system n for the non-invariant function u. From the third equation of 
(15. ip it follows that u is linear on x: 

u = k{t,y)x + U{t,y). (5.2) 

At that, function k have expression in terms of invariant functions: k = —(h t + vh y )/h. 
For the sake of convenience, one can treat this expression as an additional equation of 
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the factor-system E/H. Substitution of the representation (15.21) into the equations f l5.ll) 
and splitting with respect to x leads to the system for invariant functions: 

V t + VVy +hy — 0, 

h t + vh y + kh = 0, (5.3) 

k t + vk y + k 2 = 
and to equation for function U. 

U t + vU y + Uk = 0. (5.4) 

Equations (15. 3ft form the factor-system E/H, whereas equations ( 15. 2ft . (15.41) represent 
trivially consistent system II for the non-invariant function. The factor-system (15.31) 
itself admit some Lie group symmetries. According to lemma [3] the admissible group 
contains the subgroup with Lie algebra 
Nor L9 {Xi,X 4 } 

1— — r^: — |A 2 , A 5 , A 10 , A 11; A 12, A 13 J. 

-(Ai, A4J 

In particular, this algebra contains the subalgebra N/H = {X 10 +X 12 }. For construction 
of N/ if -invariant solution of the if -PIS, operator X10+X12 should be re- written in terms 
of invariants of algebra H in the following form (for convenience, it is also prolonged to 
the invariant variable k) 

(t 2 + l)d t + tydy + {y- tv)d v - 2thd h + (1 - 2tk)8 k . 

Invariants of this operator are 

A = yVt 2 + 1, V = vVt 2 +T-t\, H = h(t 2 + 1), K = k{t 2 + 1) - t. 

Necessary condition of an invariant solution existence is obviously satisfied. 
Representation of the invariant solution of the factor-system H 5 . 3 f) have the form: 

V(X) + tX H(X)_ K(X)+t 

v ~ Vt^+T ' t 2 + V k ~ t 2 + l ■ (5 ' 5) 

Substitution into the equation gives 
VV + H'= -A, 

VK' + K 2 + 1 = 0, (5.6) 
VH' + HV = -KH. 

Equations (15.61) form a factor-system (E/H) /(N/H) = E/N for the iV-PIS of equations 
E. The corresponding system II is given by expression (15.21) and equation (15 ,4p with 
substitution of invariant functions ( 15.51) . For the integration of the obtained system new 
independent variable /1 is introduced: 

dX t r / \ \ f dX , r _ N 

Then, the second equation of (15. 6p accurate to insufficient constant give 

K = -tgfx. (5.8) 
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By using (15. 8p . one can integrate equation (15.41) : 
g(n - aidant) 

u = TttTT (5 ' 9) 

cos /iVi + 1 

(g is an arbitrary function). Besides, system (15.61) has a first integral, which follows 
from its third equation 

HV cos \i = m. (5.10) 

Finally, system (I5.8P possess a Bernoulli integral 

V 2 + A 2 + AH = b Q . (5.11) 

The latter should be observed as an implicit (not resolved with respect to the derivative) 
equation for the dependence X(p): 

(A') 2 + A 2 + = 6 . (5.12) 

A' cos p 

Hence, the A-PIS is finally given by expressions (I5.2p . (I5.5p . where functions V, H, K, 
and U can be found form (I5.7p - (l5.1ip after integration of the first-order ODE (I5.12p . 
Note, that the solution contains an arbitrary function g. 

Analysis of the optimal system [21] for the 9-dimensional Lie algebra Lg [2J, 
admitted by equations (15.11) shows that the combination of operators {d x , td x + d u } 
or the equivalent combination {d y , td y + d v } is presented in 9 three-dimensional 
representatives. All PISs of defect 1 generated by these subalgebras are decomposable 
and can be obtained by invariant reduction of equations (15. 3p by means of one of the 
following operators: 

Xi, An, Aio + An, A5 + A10, A10, aAn + A13, A5 + An + A13, 
aAn + Aio + A12, 

Here a is an arbitrary real parameter. 



6. MHD with general state equation 

Equations of ideal magnetohydrodynamics [22l [23] are observed: 

D p + pdivu = 0, 

D u + p^Vp + p _1 H x rot H = 0, 

Dp + A(p, p) div u = 0, 

D H + H div u - (H • V)u = 0, 

divH = 0, D = d t + u-V. 

Here u = {u, v, w) is a velocity vector, H = (if, K, L) is the magnetic field; p and p 
are pressure and density. Thermodynamical functions are related by the state equation 
p = F(S,p) with entropy 5*. Function A(p, p) is determined by the state equation as 
A = p [OF /dp). All functions depend on time t and Cartesian coordinates x = (x, y, z). 
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The admissible group is 11-dimensional Galilean group extended by homothety 
[SI 121] ■ Infinitesimal operators form Lie algebra L n with basis 

X\ = d x , X 2 = d y , A 3 = d z , A 4 = td x + d u , A 5 = td y + d v , A 6 = td z + d w , 

X 7 = yd z - zd y + vd w - wd v + Kd L - Ld K , A 8 = zd x - xd z + wd u - ud w + Ld H - Hd L , 

X 9 = xd y - yd x + ud v - vd u + Hd K - Kd H , X 10 = d t , X u = td t + xd x + yd y + zd z . 

Optimal system of subalgebras QLu was constructed in [2SI [2S]; in the final form it 
can be found in [27]. Regular partially invariant solutions for gas dynamics equations 
generated by representatives of this optimal system were investigated in papers [8]-|llj. 
By virtue of special form of operators of Lie algebra L n , partially invariant solutions 
of MHD are generated by the same representatives of the optimal system as in pure 
gas dynamics without the magnetic field. Hence, the known results can be taken into 
account in construction of the solutions for MHD equations. Only non-barochronous 
solutions will be observed below, i.e. such solutions that pressure p depends on spatial 
variables x. Barochronous motions of fluid (pressure depend only on time) were carefully 
observed in general form in [T3] [15]. 

As long as all operators of Lie algebra L n have nonzero coordinates at 
differentiations with respect to independent variables, one- dimensional subalgebras do 
not give rise to a PISs. The only two-dimensional subalgebra, which generate PIS of 
defect 1 and rank 3 is {Xi,X^\. 

Analysis of three-dimensional representatives of the optimal system shows that 
there are 15 subalgebras, responsible for PISs of defect 1 and rank 2. All of 
these solutions are generated by indecomposable PISs constructed on the following 
subalgebras 

{X 2 , X 3 , X 7 }, {A 5 , Xq, X 7 }, {X 7 , Xs, X 9 }, 

(6.1) 

{A 3 + A 5 , A 2 — A 6 , A 7 }, {A 3 , A 5 , A 2 + A 6 }. 

All the rest of submodels of defect 1 and rank 2 can be obtained by the invariant 
reduction of {Ai, A 4 }-PIS with respect to one of the following operators 

X 7 + aXu, X 7 + A 10 , aA 6 + An, A 5 + A 10 , A 10 , A 2 + A 6 , A 6 , A 2 . 

There are 31 four- dimensional representatives of the optimal system, which give 
rise to PISs of defect 1 and rank 1. These are representatives of QLn with numbers 1, 
4, 5 (at a = 0), 6, 7 (at a = 0), 9 (at (5 = 0), 10 (at a = 0), 12-14, 16 (at a = 0), 17-21, 
23, 29, 30, 35, 36, 38, 41-46, 48, 49 (numeration is given according to [27]). There are 
only 9 indecomposable PISs among them with bases 

{X U X S , A 6 , «A 4 + A 7 ; a + 0}, {aX 1 + A 4 , A 5 , A 6 , (3X 1 + A 7 ; (3 + 0}, 

{X U X 2 , A 3 , «A 4 + A 7 ; a + 0}, {aX 1 + A 4 , A 3 + A 5 , A 2 - A 6 , f3X l + A 7 ; (3 + 0} 

{A 2 , A 3 , A 4 , X x + A 7 }, {A 2 , aXt + A 3 , X 1 + A 5 , A 6 ; a + 0}, 

{X\, A 3 + A 5 , A 2 - A 6 , aA 4 + I 7 ;«^ 0}, {X x , A 2 , A 3 + A 5 , A 6 }, 

{X u aX 2 + pX 3 + A 4 , <xA 3 + A 5 , rA 2 + A 6 ; a 2 + (3 2 = 1, a 2 + r V 0, P 2 + a V 0}. 
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All these subalgebras generate barochronous submodels, since their only invariant 
independent variable is time t. The only non-barochronous indecomposable partially 
invariant submodel is generated by a four-dimensional subalgebra {X 2 , A3, A5, X§}. The 
latter gives PIS of defect 2 and rank 2. 

Among five- dimensional subalgebras of Ln the only regular indecomposable and 
non-barochronous solution is generated by a subalgebra {A 2 , X 3 , A 5 , A 6 , X 7 }. For MHD 
equations this solution have defect 3 and rank 2. All the rest of regular partially invariant 
solutions generated by higher-dimensional subalgebras of Ln are either barochronous 
or decomposable. 

The calculations above show that class of indecomposable regular non-barochronous 
PIS for MHD equations is exhausted by the submodels generated by a subalgebra 
{Ai,A 4 } (defect 1, rank 3), subalgebras ( 16. II) (defect 1, rank 2), subalgebra 
{A 2 , A 3 , A 5 , A 6 } (defect 2, rank 2), and subalgebra {A 2 , A 3 , A 5 , A 6 , A 7 } (defect 3, rank 
2). Note that partially invariant solutions generated by subalgebras {A 7 , A 8 ,A 9 } and 
{A 2 , A 3 , A 7 } were already studied in [28l [29j ED] • 
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